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ABSTRACT 

Numerical simulations of outflows formed during the collapse of 100 Mq cloud cores 
are presented. We derive a generalised criterion from magnetohydrodynamical wind 
theory to analyse the launching mechanism of these outflows. The criterion is suc- 
cessfully applied to the whole outflow structure and cases with sub-Keplerian disc 
rotation. It allows us to decide whether an outflow is driven centrifugally or by the 
toroidal magnetic pressure. We show that quantities such as the magnetic field line 
inclination or the ratio of the toroidal to poloidal magnetic field alone are insufficient 
to determine the driving mechanism of outflows. By performing 12 simulations with 
variable initial rotational and magnetic energies we are able to study the influence 
of the initial conditions on the properties of outflows and jets around massive proto- 
stars in detail. Our simulations reveal a strong effect of the magnetic field strength 
on the morphology of outflows. In runs with weak fields or high rotational energies, 
well collimated, fast jets are observed whereas for strong fields poorly collimated, low- 
velocity outflows are found. We show that the occurrence of a fast jet is coupled to 
the existence of a Keplerian protostellar disc. Despite the very different morphologies 
all outflows are launched from the discs by centrifugal acceleration with the toroidal 
magnetic field increasingly contributing to the gas acceleration further away from the 
discs. The poor collimation of the outflows in runs with strong magnetic fields is a 
consequence of the weak hoop stresses. This in turn is caused by the slow build-up of 
a toroidal magnetic field due to strongly sub-Keplerian disc rotation. The mass and 
momentum outflow rates are of the order of 10~^ Mq yr~^ and 10~^ Mq km s~^ 
yr~^, respectively. The mass ejection/accretion ratios scatter around a mean of 0.3 
in accordance with observational and analytical results. Based on our results we sug- 
gest an evolutionary scenario for the earliest stage of massive star formation in which 
initially poorly collimated outflows develop which successively get better collimated 
during their evolution due to the generation of fast jets. 

Key words: hydrodynamics - magnetic fields - methods: numerical - stars: formation 
- stars: massive - stars: winds, outflows. 



1 INTRODUCTION 

The formation of massive stars and the evolution of asso- 
ciated protostellar outflows is still a highly debated ques- 
(e, 



tion 



Beuther & Shepherd 



2005 



Zinnecker & Yorke 



2007). It is believed that massive stars form in high-mass 



molecular cloud cores with masses ranging from about 100 
Mq up to a few 1000 Mq. Such cores have characteris- 
tic sizes o f a few 0.1 pc and peak densities of up to 10^ 



(e.g. Beuther et al.|2007 ). A crucial ingredient for the 
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formation of outflows is the magnetic field in the cloud cores. 
Its importance can be estimated by the mass-to-flux ratio 
fi normalised to the critical mass-to-flux ratio ([Mouschovias| 



& Spitzer 1976| which is essentially the ratio of gravita- 
tional to magnetic energy in the region (see equation ([6|). 
For high-mass star forming cores the observed mass-to-flux 
ratios are typically slightly supercritical with /x < 5 (|Fal-| 
garone et al.||2008| [Girart et al.|[2009l [Beuther et al.||2010| 



Crutcher et al. 2010) indicating a significant influence of 



magnetic fields on the star formation process. In turbulence 
simulations, however, also higher values of /i up to 20, i.e. 



weaker magnetic fields, have been found (e.g. Padoan et al 
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2001| [Tilley Pudritz||2007t . In co mbination with the ob- 



served overall slow rotation of cores (Goodman et aL||1993 



Pirogov et al.|2003t all necessary ingredients for the forma- 
tion of protostellar outflows are present. Indeed, there are a 
growing number of observations of outflows around massive 
protostellar objects (see Beuther et al.|[2002b Zhang et al 



2005 for recent compilations). The generation of such mas- 
sive outflows, in particular their underlying driving mecha- 
nism and their properties will be the focus of our interest 
here. 

For low-mass star formation (Mcore ^ 1 Mq), the gener- 
ation, evolution and properties of protostellar outflows have 
been studied in great detail over the last years (e.g. Allen 



et al.| 20031 |Banerjee & Pudritz |2006| |Mellon & Li„2008^ 


Machida et al.| 


20081 IHennebelle & Fromang| |2008| |Hen- 


nebelle & Ciardi 


2009 : Commercon et al.|2010||Tomida et al. 


2010 Duffin & Pudritz 2009; Dufiin et al. 


2011). Despite 



the intensive research in this field there is still no consensus 
about what drives the outflows beside the fact that all out- 
flows are magnetically driven. Different results can be found 
in literature about whether the outflow is driven by cen- 



trifugal accele ration (Blandford Payne 1982 Pudritz & 
Norman||T986l [Pelletier k Pudritz||1992fQr by the pressure 



of the toroidal magnetic field ( |Lyiiden-BelI||1996| |2003| ). A 

possible reason for this might be the different methods used 
to analyse the outflows. In this work we seek to clarify this 
problem by deriving a fully self-consistent, generalised crite- 
rion from magnetohydrodynamical (MHD) wind theory to 
determine the underlying driving mechanism. The criterion 
is applicable in the same way to low- and high-mass proto- 
stellar objects. 

Numerical studies on the influence of magnetic fields in 
massive star forming regions have received attention only re- 
cently. 'Banerjee & Pudritz ( 2007 ) study the formation of a 



protostar and the associated outflow in its earliest evolution 
with extremely high spatial resolution observing an outflow 
consisting of two components. The interplay of magnetically 
driven outflows and ionising radiation is analysed by [Peters] 
et al. (2011) who found a sphere-like outflow morphology. 



In this context we also mention the work of [Vaidya et al.| 
(2011) who study the effect of outflow decollimation due 



to radiation forces, but at much later times than presented 
here. Examining the influence of initial turbulence, an as- 
pect not considered in this work, Hennebelle et al. (2011) 



find that strong magnetic fields reduce the number of frag- 
ments formed during the collapse. In an accompanying pa- 
per ( Seifried et al.|20TT in the following Paper I) we analyse 
the influence of the initial condition in massive star forming 
regions on accretion rates and the build-up of protostellar 
discs. One of the main results of this study are the remark- 
able similar accretion rates despite large variations in the 
initial conditions which we attribute to effects of the mag- 
netic field. A second crucial result which will be of particu- 
lar importance for this paper is the fact that for simulations 
with strong magnetic fields (/i < 10) only sub-Keplerian pro- 
tostellar discs can form due to the very efficient magnetic 
braking. 

In the present study, we systematically analyse the in- 
fluence of the rotational and magnetic energy on protostel- 
lar outflows. We focus on the underlying launching mech- 
anism and how the initial conditions affect global outflow 
properties like mass, momentum and morphology. In par- 



ticular, we are able to examine how the development of 
a two-component outflow (Banerjee & Pudritz 2007| or a 
sphere- like outflow ( Peters et al.||201 1) depends on the ini- 
tial configuration. For this purpose we use the same simula- 
tions as in [Papei^ We note that the initial conditions are 
selected in a way to cover a large parameter space in accor- 
dance with observations and numerical simulations. With 
the generalised outflow criterion derived in this work we 
will show that magneto-centrifugally driven outflows con- 
sist of two different regimes described by the MHD wind 
theory. In the first regime close to the disc and the rotation 
axis, acceleration is dominated by the centrifugal force, i.e. 
gas gets flung outwards along the poloidal magnetic field 
lines, whereas in the second regime further away from the 
disc Bel) starts to dominate the acceleration. Furthermore we 
will also show that despite large morphological differences 
all outflows observed in this work are launched by centrifu- 
gal acceleration and that the differences in coUimation are 
mainly caused by varying hoop stresses. 

The paper is organised as follows. In Section [2] the nu- 
merical techniques and the initial conditions of the runs per- 
formed are described briefly. Next we derive the generalised 
criterion to determine the outflow driving mechanism (Sec- 
tion l3|. The results of our simulations are presented in Sec- 
tion^ First, two representative simulations are discussed 
in detail also testing the newly developed outflow criterion. 
Next, we analyse the long term evolution and stability of the 
outflows before the effect of varying initial conditions on the 
outflow properties are examined. In Section |5] the results are 
discussed in a broader context and are compared to theoret- 
ical, numerical and observational studies before we conclude 
in Section [H 



2 NUMERICAL TECHNIQUES AND INITIAL 
CONDITIONS 

In the following we briefly describe the initial conditions 
and the numerical methods used in our simulations. For a 
detailed description we refer the reader to [Paper l[ The 3- 
dimensional, MHD simulations of collapsing cloud cores are 
performed with the AMR code FLASH ( Fryxeh et al. 2000) 
using a MHD solver developed and tested by |Waagan et al.| 
( [2QTT] ). The maximum spatial resolution in the simulations 
is Ax = 4.7 AU and the Jeans length is resolved everywhere 
with at least eight grid cells. We use a threshold density of 

Pcrit = 1.78 • 10~"^^gcm~^ 



(1) 

for sink particle creation (for further details on the creation 
criteria see | Federrath et aL]|2010 l). 

With sink masses of up to Mmax 4 Mq , the maximal 
rotation velocities in the discs and outflow velocities directly 
coupled to them ( |Michel|1969| [Pelletier &; Pudritz||1992 ) 



'^out ^ '^rot,n 



Ax 



(2) 



are limited to about 10 km s ^. 

The thermodynamical evolution of the gas is modelled 



by a cooling routine of [Banerjee et al. (2006) taking into 
account dust cooling, molecular line cooling and the effects 
of optically thick gas. After the launching of the outflows, 
we introduce a density threshold of 1 • 10~^^ g cm~^ within 
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67 AU around the simulation centre. This helps to limit the 
Alfvenic velocity 

B 



(3) 



thereby preventing prohibitively small timesteps. As shown 
in IPapefT) the dynamical influence of this artificial density 
threshold is negligible. 

The simulated cloud cores have a mass of 100 Mq 
and a uniform temperature of 20 K. This mass corresponds 
to about 56 Jeans masses making the whole configuration 
highly gravitationally unstable. The core has a diameter of 
0.25 pc and is embedded in an ambient medium having a 
100 times lower density than at the edge of the core and a 
100 times higher temperature of 2000 K to assure pressure 
equilibrium. The core has a density profil^ 



p{r) oc r 



(4) 



We note that for the purely hydrodynamical case, this slope 
defines the transition between systems with low degree of 
fragmentation (for steeper density profiles) and with high 
susceptibility for fragmentation (for shallower density pro- 
files) as shown by [Girichidis et aL (2011). 

Whereas the mass distribution described above is kept 
fixed, the rotational and magnetic energies are changed be- 
tween the individual simulations. The cores are threaded by 
a magnetic field with variable strength parallel to the z-axis 
and declining outwards with the cylindrical radius R as 



oc R~ 



(5) 



Bz does not change along the 2;- axis to guarantee VB = 0. 
In addition, the cores are rotating rigidly with the rotation 
axis parallel to the initial magnetic field. We have performed 
12 simulations with varying rotational and magnetic ener- 
gies. The individual runs with the corresponding parameters 
are listed in Table ^ Additionally, the initial values of all 
runs are shown in Fig.[l] Rotational (/3rot) and magnetic en- 
ergies (7) are normalised to the gravitational energy. As can 
be seen, both /3rot and 7 cover a wide range over more than 
two orders of magnitude around the line of equipartition, 
i.e. Prot = 7- For comparative purposes the corresponding 
mass-to- flux ratio ( Mouschovias Spitzer|1976 ) 



M ■ 



'/(-) 



0.13 



(6) 



and the angular frequency u are listed. We note that the 
inclusion of turbulent motions in the setup will be delayed 
to a subsequent paper. 



3 A GENERALISED CRITERION FOR 

MAGNETO-CENTRIFUGALLY DRIVEN 
JETS 

Several approaches can be found in literature to deter- 
mine whether an outflow is driven by centrifugal acceler- 
ation (iBlandford & Paynell 19821 iPudritz k Norman||1986 



^ To avoid unphysically high densities in the interior of the core, 
we cut off the r"-*^- ^-profile at a radius of 0.0125 pc. Within this 
radius the density distribution follows a parabola p(r) oc (1 — 



Table 1. Performed simulations with normalised initial mass-to- 
flux-ratio /i, ratio of magnetic to gravitational energy 7, magnetic 
field strength in the centre B, ratio of rotational to gravitational 
energy /3rot, and the corresponding angular frequency lj. 



Run 






7 








rin-lS -11 


26-20 


26 


1.13 


• 10-3 


132 


2 


•10-1 


7.07 


26-4 


26 


1.13 


• 10-3 


132 


4 


• 10-2 


3.16 


26-0.4 


26 


1.13 


• 10-3 


132 


4 


• 10-3 


1.00 


26-0.04 


26 


1.13 


• 10-3 


132 


4 


• 10-4 


0.316 


10-20 


10.4 


7.06 


• 10-3 


330 


2 


•10-1 


7.07 


10-4 


10.4 


7.06 


• 10-3 


330 


4 


• 10-2 


3.16 


10-0.4 


10.4 


7.06 


• 10-3 


330 


4 


• 10-3 


1.00 


5.2-20 


5.2 


2.82 


•10-2 


659 


2 


•10-1 


7.07 


5.2-4 


5.2 


2.82 


•10-2 


659 


4 


• 10-2 


3.16 


5.2-0.4 


5.2 


2.82 


•10-2 


659 


4 


• 10-3 


1.00 


2.6-20 


2.6 


1.13 


•10-1 


1318 


2 


•10-1 


7.07 


2.6-4 


2.6 


1.13 


•10-1 


1318 


4 


•10-2 


3.16 




Figure 1. Initial values of the 12 simulations performed. The 
initial rotational (/3rot) and magnetic energy (7) are normalised 
to the gravitational energy. As can be seen, the starting points 
bracket the curve where rotational and magnetic energy are equal 
(black line). The second axes show the corresponding normalised 
mass-to- flux ratio ji and the rotation frequency uj. 



Pelletier & Pudritz 1992) or the magnetic pressure gradi- 
ent (Lynden-Bell 1996 2003). A frequently chosen quan- 



tity is the ratio of the toroidal to poloidal magnetic field, 
B(f)/B^o\ (e.g. |Hennebelle & Fromang 2008). If this ratio 
is significantly above 1, the outflow is often believed to be 
driven by the magnetic pressure. However, the considera- 
tion of B(f) I Bpoi alone can be misleading as in centrifugally 
driven flows this value can be as high as 10 (see Fig. 4 in 
Blandford Payne] [19821 henceforth BP82). Close to the 
disc surface one can check the inclination of the magnetic 
field lines w.r.t. the vertical axis. The field lines have to 
be inclined by more than 30° for centrifugal acceleration to 
work (BP82V Although this criterion is an exact solution of 
the ideal, stationary and axisymmetric MHD equations for 
an outflow from a Keplerian disc, its applicability is limited 
to the surface of the disc. A criterion to determine the driv- 



ing mechanism above the disc was used by Tomisaka ( 2002 ) 



comparing of the centrifugal force Fc and the magnetic force 
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Fmag. By projecting both forces on the poloidal magnetic 
field fines it can be determined which force dominates the 
acceleration. For the outflow to be driven centrifugally, Fc 
has to be larger than Fmag. However, for this criterion to be 
self-consistent the gravitational force and the fact that any 
toroidal magnetic field would reduce the effect of Fc have to 
be taken into account. 

In the following we derive a generalised criterion based 
on the ideal, axisymmetric, time-independent (stationary) 
MHD equations ( |Chandrasekhar||1956] |Mestel|p^6Tt to de- 
cide whether centrifugal acceleration dominates or not. The 
criterion is applicable everywhere in the outfiow and not only 
at the disc surface. To begin with we review the basic results 
of the solution of the stationary, axisymmetric MHD equa- 
tions. As the solution is discussed broadly in literature (e.g. 
Blandford Payne|[T982l [Pelletier Pudritz| | 19921 |Spruit| 
1996) we only state the main results necessary to derive our 



Vcf> 

r 



(8) 



criterion. One key aspect of the solution is that it invokes 
four surface functions which are constant on each magnetic 
surface, i.e. along each magnetic field line. We list them in 
the following. 

1. Combining the induction equation, the divergence-free 
equation and the continuity equation shows that the poloidal 
velocity is always parallel to the poloidal magnetic field 

pvpoi = kBpol (7) 

where the so-called mass loading, is the first surface func- 
tion, i.e. it is constant along any given field line. 

2. The second surface function which can be derived from 
the three aforementioned equations is the angular velocity of 
the magnetic surface 

kB^ 
rp 

Here r is the cylindrical radius and Vcj) the local rotation 
speed of the gas. 

3. Making use of the equation of motion yields the third 
surface function, the angular momentum invariant 

Ank 

4. Finally, considering the kinetic energy equation yields 
the energy invariant along a poloidal field line 



2 Auk 

^2 |l2,^,7 Vet) ^ Bcf^Bpol . 1 

where h is the enthalpy and $ gravitational potential. This 
equation is also called the Bernoulli equation. 



^ rvct> ■ 



(9) 



In the second line of equation ( 10 ) we have used the 



variables which can be inferred from simulations to describe 
e. We especially note the last term containing B^ describing 
the infiuence of the toroidal magnetic field on the dynamics 
of outfiows. It states that in magnetic wind theory besides 
the poloidal magnetic field fiinging material outwards, the 
toroidal magnetic field can contribute to the acceleration of 
gas (see also |Spruit|1996| for a detail discussion). Magneto- 
centrifugal acceleration therefore has two regimes; a centrifu- 
gally dominated and a (magnetic pressure) dominated 
regime. In the following we will use centrifugal acceleration 



when describing regions in which gas is accelerated mainly 
by fiinging gas outwards along the poloidal field lines. In 
contrast, magneto -centrifugal acceleration also contains the 
possibility of dominated acceleration. 

In the following we always assume that the thermal 
pressure can be neglected, i.e. /i = 0. Whether gas can get ac- 
celerated, i.e. whether Vpoi increases, depends on the change 



of the second to last term on the r.h.s. of equation ( 10) (sec- 



ond line) along the poloidal field line where e is constant. 
The general condition for acceleration is therefore 



* 2 ^ Vpoi 47r 



1 B(i)Bpo\ _^ 1 ^0 

p 4:7T p 



< 



(11) 



where ^poi denotes the derivative along the poloidal mag- 
netic field. This criterion is universal and can be applied to 
the whole outfiow structure. It describes all regions of gas 
acceleration including those dominated by the effect of B^. 
Furthermore it only contains variables which are directly 
accessible from the simulations. 

Now, we come back to our original aim to work out a 
criterion for identifying the regions dominated by centrifu- 
gal acceleration. To do so, we find it useful to go to a frame 
rotating with the magnetic field, i.e. with the angular fre- 
quency (jj. For this purpose the Bernoulli equation (10) is 
modified as follows 



lu ■■ 



where $cg is the centrifugal-gravitational potential 



(12) 



(13) 



with M being the mass of the central object and G the grav- 



itational constant. Compared to equation (10) all magnetic 
field terms have vanished. This is due to the fact that in 
the rotating frame the gas fiow is strictly parallel to the 
magnetic field and therefore the Lorentz force vanishes. 

We now make the additional assumption that the mag- 
netic field is strong enough to retain a purely poloidal struc- 
ture, i.e. ^0 = 0, and to force the gas to corotate with the 
field. Hence, we can set v^/r = cj (compare equation (|8|) 
causing the second term on the r.h.s. of equation (12) to 
vanish. This in turn reduces the question of whether gas 
can be accelerated to the question of how $cg changes along 
the field line. $cg decreasing along the field line causes I'poi 
to increase, i.e. the gas gets accelerated and vice versa. We 
emphasise that this works because gas can only move paral- 
lel to the poloidal magnetic field in the stationary case (see 
equation ([9|). We note that setting uj in equation ( 13 ) to the 
Keplerian value at the footpoint of the magnetic field line 
in the disc yields the well-known criterion found by ,BP82| 
stating that the field lines at the disc surface have to be 
inclined by more than 30° w.r.t. the z-axis for centrifugal 
acceleration to work. However, the expression here is more 
general and sub-Keplerian velocities may also be used. 

In a numerical simulations, for an arbitrarily chosen 
point somewhere above the disc it is often not possible 
to determine the footpoint of the magnetic field line pass- 
ing through that point. Therefore, it is not possible to set 
UJ = where is the angular frequency at footpoint in 
the disc where the field line is anchored. Under the assump- 
tion that gas and magnetic field corotate, we can circumvent 
this problem. We therefore replace uj by the local angular 
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frequency of the gas, i.e. cu — Q = t'(/,/r which is known 
from the simulation data. Hence, at any given point (R*,Z*) 
we can calculate the value of the centrifugal-gravitational 
potential 



of corotation, i.e. = 0, made for the derivation of equa- 



tes,* $cg(R*,Z*) 



CM 



+ Z2 2 



1 2 



(14) 



By solving the equation 



ixj^r — <Pcg,* = const. [Id) 



+ ^2 2 

under the assumption that the gas would rotate everywhere 
with cj* = one obtains the shape of an isocontour 



z(r, cj*, $cg,*) 



r2cj2 + 2$c 



(16) 



Along this line the value of $cg stays constant, i.e. $cg = 
$cg,*- Assume the magnetic field line passing through the 
point (R*,Z*) is rigidly corotating with the gas with the an- 
gular velocity cj* , one can easily decide whether the gas can 
get accelerated or not by comparing the shape of the iso- 
contour with that of the magnetic field line. Hence, taking 
the derivative of equation (16) w.r.t. r and comparing it to 
the ratio Bz/Br at the point (R*,Z*) yields the desired cri- 
terion. For centrifugal acceleration to work at an arbitrarily 
chosen point the criterion therefore is 



r 1 

z GM 



-^0 / 2 I 2n3/2 



GM 



(I) 



> 1 



(17) 



For the sake of simplicity we have neglected the * symbol in 
the above equation. Every quantity has to be taken at the 
point considered, i.e. at (r, z). 



We emphasise that the criterion given in equation (17) 



traces the regions where acceleration is mainly by the cen- 
trifugal force, i.e. flinging out gas along the poloidal field 
lines. As we have set = 0, there is no resulting Lorentz 
force along the poloidal field line (see also Section 4.1 of 
Spruit p 996) and the only acceleration is indeed due to the 
centrifugal force exceeding the gravitational force. The mag- 
netic field only enforces the gas to rotate with the angular 
frequency cj. Regions of acceleration not traced by this cri- 



terion should be fit by the criterion given in equation (11) 
as here the acceleration by is taken into account. There- 
fore, by comparing the results of both criteria we are able 
to distinguish between regions dominated by centrifugal ac- 
celeration and those by the toroidal magnetic pressure. 

We again emphasise that the formulation given in equa- 



tion (11) states that in magnetic wind theory, the acceler- 
ation by B(f) in magneto-centrifugal flows is implicitly con- 
tained. This seems often to be overlooked in literature when 
discussing the underlying outflow mechanism. Equation [TT] 
nicely demonstrates that a classification into centrifugally 
driven flows on the one hand and magnetic tower flows on 
the other might be an oversimplification as in reality there 
can be a coexistence of centrifugal acceleration and acceler- 
ation due to the pressure gradient of B^. We will refer to 
this several times in the paper. We also point out that the 
equations [To] and [12] are equivalent despite their very dif- 
ferent appearance. The equations and |17| however, are 
not equivalent any more due to the additional assumption 



tion (17) 



4 RESULTS 

In this section, we present the results of our collapse simu- 
lations focussing on the evolution and the launching mecha- 
nism of outflows. The evolution of protostellar discs and the 
accretion properties of the protostars were studied in detail 
injPapeFT] We limit our consideration to the phase after the 
first sink particle has formed to study the time evolution 
of the outflow and its underlying launching mechanism in 
detail. In the following, all times refer to the time elapsed 
since the formation of the first sink particle. 

In general, the outflows are launched shortly after the 
creation of the first sink particle as soon as a protostellar disc 
builds up. In each simulation the outflow evolves over time 
and shows no signs of re-collapse until the end of the simu- 
lation. However, outflow morphologies and global properties 
like mass or momentum (see Table [5]) differ significantly be- 
tween the individual runs. Hence, in the following we select 
two representative simulations with equal initial rotational 
energies but different magnetic field strengths to study the 
properties and the launching mechanism of the outflows in 
detail. 



4.1 Weak field case 26-4 

4.1.1 General properties 

We start our consideration with run 26-4 which has a weak 
initial magnetic field and a moderate rotational energy (see 
Table[T]). In Fig.[2]we show the density, the poloidal velocity 
field and the outflow velocity 



\z\ 



(18) 



of the outflow in a slice along the z-axis at two different 
times. It can be seen that initially the outflow expands rather 
slowly into the surrounding medium as within the first 2000 
yr it has reached a height of only ^ 250 AU. After 2000 
yr the growth rate increases and reaches a roughly constant 
value of about 1 AU yr~^ ^ 4.7 km s~^. We will discuss 
this behaviour in detail at the end of Section |4.1.2| Fig. |2] 
also shows that a bow shock at the tip of the outflow de- 
velops extending down to the edge of the centrifugally sup- 
ported disc. As indicated by the velocity vectors, most of 
the gas within the bow shock is directed vertically outwards 
with velocities up to 15 km s~^. This corresponds to highly 
supersonic motions of Mach numbers up to ~ 55 using a 
temperature of 20 K as present in the undisturbed, ambient 
medium and Mach numbers of about 15 - 20 w.r.t. tem- 
peratures of a few 100 K in the outflow itself. The outflow 
morphology at 2000 yr and 5000 yr reveals a self-similar ap- 
pearance. In particular the collimation factor of the outflow, 
i.e. the ratio of the length to the width of one outflow lobe, 
settles around a value of 4. We especially note the very 
turbulent structure seen in the outflow after 5000 yr in both 
the density field and the outflow velocity. Several internal 
shock fronts and instabilities have occurred in the outflow 
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Table 2. Listed are the mass, momentum and extension of the outflows, the total sink particle mass at the end of each simulation, the 
time the runs have being followed after the flrst sink particle has formed, and the time averaged mass and momentum outflow rates. 



Run 


Mout 


^out 


L 


Msink 


tsim 


Mout 


^out 




(Mq) 


(M0 km s-i) 


(AU) 


(M0) 


(yr) 


(10-4 M0 yr-i) 


(10-4 M0 km s-i yr-i) 


26-20 


0.509 


0.337 


707 


1.85 


5000 


1.02 


0.673 


26-4 


0.853 


2.37 


3215 


2.65 


5000 


1.71 


4.74 


26-0.4 


0.526 


2.23 


3720 


3.59 


5000 


1.05 


4.46 


26-0.04 


0.080 


0.202 


856 


4.26 


5000 


0.16 


0.404 


10-20 


1.09 


1.02 


1240 


1.28 


4000 


2.73 


2.54 


10-4 


0.603 


0.585 


669 


2.23 


4000 


1.51 


1.46 


10-0.4 


0.164 


0.177 


445 


2.98 


4000 


0.41 


0.442 


5.2-20 


0.875 


1.70 


2110 


1.78 


4000 


2.19 


4.25 


5.2-4 


0.656 


0.715 


1128 


2.28 


4000 


1.64 


1.79 


5.2-0.4 


0.537 


0.586 


1100 


2.55 


4000 


1.34 


1.46 


2.6-20 


0.116 


0.281 


1942 


1.30 


3000 


0.39 


0.938 


2.6-4 


0.095 


0.215 


1455 


1.48 


3000 


0.32 


0.717 



although no perturbations are included in the initial condi- 
tions. This turbulent structure is a consequence of instabil- 
ities and not of episodic ejection events. We checked this by 
visual inspecting the time evolution of the outflow, finding 
that continuous ejection occurs over the whole time range. 
Furthermore, a highly time variable ejection rate would also 
require significant variations in the mass accretion rate (e.g. 
Pudritz Normaii||1986| ) which is clearly not the case (see 
Fig. 13 in Paper I). A turbulent structure as observed here 



was recently also reported by Staff et al. (2010) for jet simu- 



lations. The energy in the shocks dissipates and heats up the 



jet possibly resulting in optical forbidden line emission (Staff 
et al.|2010| . Interestingly, despite the turbulent structure the 



outflow keeps expanding with constant speed. Another in- 
teresting fact is that the gas is continuously ejected from the 
disc over the whole 5000 yr although the protostellar disc 
starts to fragment around t ^ 2500 yr and has formed 12 
further fragments by then (see Paper I[ for details). Hence, 
fragmentation does not necessarily terminate the driving of 
outflows, a fact which is also observed in run 26-20. 

To study the outflow structure in more detail, we show 
the position- velocity (PV) diagram of the outflow at 5000 yr 
m Fig. |3] The PV diagram is frequently used by observers 



to study the velocity structure of outflows (e.g Lada & Fich 
1996| [Beltran et aL]|2011t . For the dia gram shown here we 



only take into account outflowing gas. As indicated by the 
contours, the velocity of the bulk of outflowing material in- 
creases with distance from the midplane. Such a "Hubble 
Law" is frequently observed for outflows around low- and 
also high-mass protostell ar objects (e.g. p^ada &: Fich||1996| 
Arce Goodman||2001| [Beuther et al.(|2003| |Wang et ah] 



2011||Ren et al.||2Qll| [Beltran et aL||2Qllt . The maximum 



outflow velocity saturates around 15 km s~ showing several 
distinct peaks which we attribute to the internal shocks in 
the outflow (see bottom panel of Fig.|2]). Above a distance of 
about 2000 AU, however, the maximum velocity experiences 
a significant drop from about 15 to 10 km s~^. At the same 
time the "Hubble law" describing the evolution of the bulk 
velocity truncates and the bulk velocity saturates at a value 
of ^ 5 km s~^. We attribute this to a strong internal shock 
front and the very turbulent structure at these distances 



which reduce the maximum outfiow velocity and prevent an 
efficient overall acceleration although local gas acceleration 
is still possible. The position where this happens gradually 
moves outwards whereas the bulk velocity above this point 
remains almost constant with a value of ^ 5 km s~^. 

In Fig. [4] we show radial profiles of the density and out- 
fiow velocity at different vertical positions. Both quantities 
are averaged azimuthally before plotting. As can be seen, 
below 2000 AU the highest velocities occur close to the sym- 
metry axis of the jet. In particular the velocity profile at z 

— 500 AU shows a very well confined, fast velocity compo- 
nent with a quite sharp drop-off at a radius of 150 - 200 AU. 
This two velocity components, the fast, central as well as 
the slower, enclosing component are also seen in the PV di- 
agram (Fig. [3| in particular in the region between z — and 

— 1000 AU where two "Hubble laws" seem to be present. The 
good coUimation of the fast jet component is most likely due 
to the strong hoop stress exerted by the toroidal magnetic 
field. Such two velocity component outfiows with a well col- 
limated, fast component close to the axis of the outfiow and 
a wider spread, slower component are frequently observed 
around low-mass protostars (see e.g. Bachiller||1996, for a 
review) and also around massive protostars (e.g. Beuther 



et al. 2004 Ren et al. 2011). Furthermore, they are also 



observed in jet simulations (e.g. Staff et aL]|2Q10| 



We note that the decrease of Vz,out close to r = at z = 
500 AU (black line in the top panel of Fig. [4| is most likely 
a numerical issue as gas very close to the z-axis cannot get 
accelerated properly due to the limited resolution. Above z 
— 2000 AU the velocity profile is more or less smooth over 
the whole radial range. In accordance with the PV diagram 
(Fig. [3} for those distances no velocities higher than ^ 6 
km s~ occur. The density profiles (bottom panel of Fig. ^ 
show a relatively fiat shape at all distances with variations 
of about half an order of magnitude and possibly a slight 
increase towards the z-axis. At the outer edge the density 
experiences a strong increase due to the bow shock confining 
the outfiow. By carefully comparing with the top panel of 
Fig. ^ it can be seen that the material in the bow shock 
partly reveals infall velocities. 
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Figure 2. Slice along the z-axis in the weak-field run 26-4, at 2000 
yr and 5000 yr after the formation of the first sink particle. The 
two top panels show the density field and the poloidal velocity 
vectors (black arrows). Note the different spatial scales between 
the top and middle panel. The outfiow velocity (^'z,out, bottom 
panel) and the density field after 5000 yr show a very turbulent 
structure caused by internal shocks. 




-3000 -2000 -1000 1000 2000 3000 
z/AU 

Figure 3. Position-velocity diagram for the weak-field run 26-4 
after 5000 yr. The contours have a logarithmic spacing. The bulk 
velocity increasing with distance. The maximum velocity show 
several clear peaks which are attributed to internal shock fronts. 
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r/AU 

Figure 4. Radial profiles of the outfiow velocity (top) and the 
density (bottom) at different vertical positions for the weak-field 
run 26-4 after 5000 yr. The quantities are averaged azimuthally 
before plotting. For z < 2000 AU the outfiow velocity increases 
towards the 2;-axis whereas at larger distances it has an almost 
fiat radial profile. The density profiles are rather fiat showing only 
small variations and a prominent jump associated with the bow 
shock. 
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4.1-2 Launching mechanism 

Next, we analyse the underlying launching mechanism of the 
outflow in the weak field case 26-4. Once again, in order to 
smooth out local perturbations which inevitably would oc- 
cur in an arbitrarily chosen slice along the ^-axis, the quanti- 
ties shown in the following are averaged azimuthally. Firstly, 
the relative importance of the toroidal (5^) and poloidal 
magnetic field (i^poi) for the total magnetic energy content 
in the outfiow is studied. This is done in the top panel of 
Fig. [5] where the value of B^/Bpoi in a slice along the z-axis 
is shown. As can be seen, almost the complete outfiow re- 
gion is dominated by i?^ . The toroidal component is created 
by the rotation of the protostellar disc in which the initially 
purely poloidal magnetic field is anchored. 5^ reaches values 
larger than Bpoi by up to 2 orders of magnitude. Hence, i?^ 
should have a crucial effect on the evolution of the outfiow. 
Depending on the position, the absolute values of i?^ and 
Bpoi vary between ^ 0.01 G and ^ 1 G. 

Outfiows driven by the pressure gradient of the toroidal 
magnetic field are described by "Lynden-Bell| ( 1996 2003) 
and are often termed as magnetic tower fiows. The fact that 
the outfiow region is mainly dominated by suggest that 
the outfiow is such a magnetic tower fiow. However, a closer 
inspection reveals that the region close to the symmetry 
axis, where the highest velocities occur (see bottom panel 
of Fig. [2]), is either only weakly dominated by the toroidal 
magnetic field or even dominated by Bpoi, i-e. B^/Bpoi ^ 10. 
We emphasise that even a field configuration which is dom- 
inated by B(f) does not contradict the centrifugal accelera- 
tion mechanism. |BP82 find values of B^/Bpoi up to ^ 10 
(see their Fig. 4) in the acceleration region in agreement 
with our findings. Therefore, we assume that the outfiow is 
driven by centrifugal acceleration although there is certainly 
a contribution to the acceleration by as well. Of course, 
this assumption has to be confirmed quantitatively in the 
following. 

To do so, in the bottom panel of Fig. [5] we analyse the 
magnetic field structure in detail concentrating on the in- 
ner region where the jet is launched. The poloidal magnetic 
field lines are overplotted on the density field. In addition, 
we show the poloidal velocity field and the contours where 
B(j)/Bpo\ = 1, 10 respectively. As can be seen, the field lines 
just above/below the disc are strongly inclined against the 
z-axis. Except for the innermost part the inclination angle is 
everywhere above 30° which is required to launch cold gas 
from the disc by centrifugal acceleration (BP82). Here we 
emphasise that an inclination angle smaller than 30° does 
not necessarily mean that centrifugal launching is impossi- 
ble. The 30° condition is valid only for a cold gas, i.e. if 
thermal pressure can be neglected. As soon as a thermal 
pressure gradient is present aiding to accelerate the gas up- 
wards, even inclination angles below 30° are sufficient for 
jet launching. As there is clearly a pressure increase towards 
the centre in the simulation, the innermost region is suit- 
able for centrifugal acceleration as well even though the in- 
clination is below 30°. Furthermore, as mentioned before in 
the inner region where the actual acceleration takes place, 
the magnetic field is not or only weakly dominated by 
in agreement with th e findings of |BP82t In contrast, for a 
magnetic tower fiow (| Lynden-Bell||2003| ) one would expect 
a much more wound up structure with a clearly dominating 
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Figure 5. Magnetic field lines structure of the weak-field run 
26-4, after 5000 yr in a slice along the 2;-axis. All quantities are 
averaged azimuthally. Top: Ratio of the toroidal to poloidal mag- 
netic field. Black contours give the transition from the toroidally 
to the poloidally dominated region. Almost the whole outflow area 
is dominated by Bottom: Density fleld in the inner region. 
Also shown are the poloidal fleld lines (white lines), the veloc- 
ity fleld (black arrows) and the contours where B^/Spoi = 1, 10 
(black and dark blue lines, respectively). Almost everywhere the 
poloidal magnetic fleld has inclination larger than 30° w.r.t. the 
z-axis suitable for centrifugal acceleration. 



toroidal magnetic field component. However, we note that 
the toroidal magnetic pressure also contributes to the accel- 
eration of the gas as implied in the MHD wind theory (see 



equation ( 11 )). As argued before, the launching from the disc 



itself, however, is most likely due to centrifugal acceleration. 

To further support the conclusion of a centrifugally 
driven wind, we apply the two criteria derive d in Section [3 



see equations 11 and 17). In contrast to the Blandford & 



|Payne| criterion, we can in general determine the driving 
mechanism of the outfiow away from the disc surface. The 
results are shown in Fig. [6] where the regions where the cri- 
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Figure 6. Slice along the z-axis in the weak-field run 26-4 after 5000 yr for two different scales. Gray shaded areas show the regions 
where the criteria derived in Section [s] (equation ( |17| (left panel) and equation (right panel)) are fulfilled. The left panels show that 
centrifugal acceleration works mainly close to the z-axis up to a height of about 800 AU which agrees very well with the region where 
the highest velocities are found. The general criterion is more volume filling and traces also regions in the outer parts. 



teria are fulfilled are shaded grey. As can be seen, within 
about 800 AU above/below the disc there is a continuous 
region close to the z-axis where centrifugally dominated ac- 
celeration is possible (left panel). Also at larger heights cen- 
trifugal acceleration partly works, in particular close to the 
symmetry axis of the jet. The grey shaded region around the 
2;-axis in the upper left panel of Fig. [6] agrees very well with 
the region where the highest outflow velocities are detected. 
Hence, despite the strong assumptions like stationarity, axis- 
symmetry and corotation used to derive the criterion given 



in equation (17), it works reasonably well. In particular, it 



is applicable to regions far from the disc surface in contrast 
to the 30° - criterion of lBP82l We therefore draw the follow- 
ing conclusions: The assumption of corotation is reasonably 
well fulfilled in the inner part of the outflow {\z\ ^ 800 AU) 
despite the partly significant toroidal magnetic field (see 
Fig. [5]). Hence, the central jet can be considered as centrifu- 
gally driven, i.e. gas gets flung out- and upwards along the 
poloidal magnetic field lines. The result also demonstrates 
that the 30° condition - apart from the disc surface - and 
the ratio of to Bpoi are not sufficient to determine the 
driving mechanism above/below the disc. 



Taking into account the effect of in the most general 
expression of magneto-centrifugal driving (equation |TT]) ) 
gives a markedly different result (see right panel of Fig. |6f . 
The grey shaded areas are now much more volume filling, 
in particular in the outer parts of the outfiow with radii 
> 200 AU where outfiowing gas is present as well. As purely 
centrifugal dominated acceleration does not work in this re- 
gion, we expect to be mainly responsible for the out- 
fiow driving. The results in the upper panel of Fig. [6] nicely 
demonstrate the capability of our criterion to distinguish be- 
tween the different driving mechanisms which was the rea- 
son for its derivation. An increasing importance of for 
the driving can also be seen in the upper parts of the outfiow 
{\z\ ^ 800 AU). Here the situation is less suitable for pure 
centrifugal acceleration (bottom left panel of Fig. [g]) but in 
general magneto-centrifugal acceleration is still possible in 
great parts (bottom right panel) . This is why we argue that 
Bcj) must contribute significantly to the outfiow dynamics at 
great heights. 

We note that even the general criterion given in equa- 



tion ( 11 ) does not cover the entire outfiow structure. This is 



due to the fact that the gas indeed does not get accelerated 
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everywhere, e.g. in the shock regions, but possibly also due 
to the azimuthal averaging process. However, the very dif- 
ferent results given by the two criteria, in particular in the 
upper panel of Fig|6] strongly indicate that the two distinct 
outflow regions are real. 

Although the gas in the outflow experiences several in- 
ternal shocks (see Fig. [2]), the bulk velocity steadily increases 
within about 2000 AU (Fig. [sf. This can be explained by 
the fact that acceleration is possible over almost the whole 
outflow extension (bottom right panel of Fig.[6|. After expe- 
riencing a shock which decreases the outflow speed, the gas 
gets reaccelerated again by centrifugal acceleration and/or 
the toroidal magnetic fleld pressure. This happens repeat- 
edly over the whole outflow extension thereby successively 
increasing the bulk velocity. This situation markedly differs 
from episodic jet ejection which also would produce inter- 
nal shock fronts. Therefore, we conclude that the knotty 
structure often observed in protostellar jets is not necessar- 
ily a consequence of several outflow ejection events but can 
also result from a continuously fed jet where gas repeatedly 
shocks and reaccelerates (see also Staff et al.|2010| ). 

As already seen in the bottom panel of Fig. [5] the mag- 
netic field lines get straightened very quickly above and be- 
low the disc. We attribute this to the hoop stress produced 
by the toroidal magnetic field collimating the outfiowing gas 
and therefore also the magnetic field lines. Nevertheless, the 
gas still gets accelerated centrifugally despite the almost ver- 
tical direction of the magnetic field lines (Fig. [gJ. Further- 
more, as shown in Fig. |5] the largest part of the outfiow is 
dominated by the toroidal magnetic field component. There- 
fore it is not surprising that over its complete extension the 
outfiow as a whole stays well coUimated. 

As mentioned in the previous section, the expansion 
speed is not constant over time but experiences a relatively 
sharp increase after about 2000 yr. We attribute this to a 
change in the underlying driving mechanism of the outfiow. 
Indeed, analysing the outfiow with our criterion given in 



equation (17) shows that within the first 2000 yr purely 



centrifugal acceleration is not possible indicating that the 
outfiow is mainly driven by the toroidal magnetic pressure. 
The expansion speed in this phase is almost the same in the 
vertical and horizontal direction (see top panel of Fig. [2| 
with the outer edge of the bubble coinciding with the po- 
sition of the accretion shock at the disc edge. It is only in 
this initial stage when we call the outfiow a magnetic tower 
fiow ( [Lynden-BeU|p^996 2003). In contrast to the situation 
at 5000 yr in this transient phase there is no acceleration 
of gas from the disc. In fact, the gas is accelerated only at 
the tip of the outfiow. Therefore the situation differs signif- 
icantly from the later stages. After ^ 2000 yr a fast, well 
collimated outfiow component, the centrifugal driven jet de- 
velops in the region close to the ^-axis. The launching of the 
jet coincides with the build-up of a well defined, extended (^ 
100 AU) Keplerian disc whereas prior to that disc rotation 
is mostly sub-Keplerian. 

In summary, besides the magnetic field line structure, 
the application of the criterion derived in this work strongly 
indicates that the outfiow is mainly driven centrifugally at 
1^1 ^ 800 AU while the dynamics of gets more important 
at large radii and larger heights where the fiow is magneto- 
cent r if ugally driven. 



4.2 Strong field case 5.2-4 

4.2.1 General properties 

Next, we describe global properties of the outfiow gener- 
ated in run 5.2-4 which has a 5 times stronger initial mag- 
netic field than run 26-4 (see Table [T]). The outfiow shown 
in Fig. [7| reveals significant differences compared to the out- 
fiow in run 26-4 (compare Fig.[2|. Whereas the latter is well 
collimated with a coUimation factor of ^ 4, the former has a 
rather sphere-like morphology expanding with roughly the 
same speed in all directions, therefore also maintaining a 
self-similar morphology for all times. The outfiow velocities 
reach values of up to 5 km s~^ about a factor of 2 - 3 lower 
than in run 26-4. The expansion speed of the outfiow is al- 
most constant over time with a value of 0.28 AU yr~^ ^ 1.3 
km s~^. This is noticeably smaller than the expansion speed 
of 1 AU yr~^ observed in run 26-4. Consequently also the 
bow shock structure in run 5.2-4 is less pronounced. 

Furthermore, a closer inspection of the outfiow in run 
5.2-4 reveals that in particular close to the symmetry axis 
and the centre of the bubble gas is still falling inwards even 
at late times. Gas with outwards motion occurs mainly in 
the outer wings. The outfiow direction in the inner part is 
almost radial and gets collimated at relatively large radii 
of ^ 500 AU. This is remarkably different to the situation 
in run 26-4 where almost all the gas within the outfiow 
area is moving outwards and preferentially parallel to the 
z-axis. A consequence of the complicated velocity structure 
observed in Fig. [7| is the complex density structure show- 
ing several shock-like features in the bubble. We find that 
the fiattened structure in the midplane is a strongly sub- 



Keplerian disc with significant infall motions (see [Paper I 
for a detailed discussion). The sub-Keplerian rotation is a 
consequence of the strong initial magnetic field decelerat- 
ing the rotation of the initial core via the magnetic braking 
mechanism ( Mouschovias Paleologou||1980 ). 



4.2.2 Launching mechanism 

The different morphologies of the outfiows in run 5.2-4 and 
run 26-4 raise the question whether the underlying launching 
mechanisms differ. Firstly, we examine the relative impor- 
tance of the toroidal and poloidal magnetic field components 
in the top panel Fig.js] Here again - as in Section [4T2I - afl 
quantities like velocity, magnetic field and density are aver- 
aged azimuthally before being plotted in order to smooth out 
local variations which might complicate the analysis. As can 
be seen, a great part of the outfiow bubble is not dominated 
by the in contrast to the outfiow in run 26-4 (see Fig.|5|. 
Interestingly, the regions where dominates are usually as- 
sociated with relatively fast outfiowing gas. In contrast, the 
regions which are dominated by Bpoi either have only slow 
outfiow velocities (< 1 km s~^) or even show infall motions. 
The fact that outfiowing gas seems to be associated with 
a strong toroidal magnetic field, i.e. B^/Bpoi > 1, suggests 
that the outfiow might be driven by the pressure gradient 
of B^. 

In the bottom panel of Fig. [S] the magnetic field line 
structure in the inner region is considered in more de- 
tail. The poloidal field lines have a strongly pronounced 
hourglass-shaped configuration with a strong radial com- 
ponent. This is caused by the rapid gas infall in the disc 
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Figure 7. Same as in Fig.|2]but for the strong-field run 5.2-4. The 
outflow is poorly collimated and has significantly lower outfiow 
velocities than the outflow in run 26-4. 



which continuously drags the field lines inwards. Although 
the poloidal field line structure seems suitable for centrifu- 
gal acceleration, i.e. the magnetic field lines are inclined by 
more than 30° w.r.t. the z-axis, only in parts of the region 
outflowing gas is present, in particular in regions more than 
200 AU above/below the disc. Hence, the 30° criterion does 
not work in this case. 
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Figure 8. Same as in Fig. [s] but for the strong-fleld run 5.2- 
4. Top: Additionally the poloidal velocity fleld (black vectors) 
is shown. Only parts of the outflow bubble, mainly associated 
with outflowing gas, are dominated by Bottom: The strong 
inclination of the poloidal magnetic fleld lines is caused by the 
inwards drag of infalling gas in the disc. 



The failure of the [Blandford Payne| criterion is not 
surprising as it was derived for Keplerian disc rotation only 
and therefore does not apply for strongly sub-Keplerian discs 
as it is the case here. The physical reason for centrifugal ac- 
celeration to fail in the inner region are the slow rotation 
velocities. As therefore the centrifugal force Fc is reduced 
significantly, close to the protostar and the disc Fc cannot 
overcome gravity and consequently gas cannot get flung out- 
wards along the poloidal magnetic field lines. Only at larger 
radii, where the gravitational force is sufficiently reduced, 
does centrifugal acceleration work. This fits with the obser- 
vation that the bulk of outflowing material indeed emerges 
from radii > 300 AU (see FigM. 



Our centrifugal launching criterion (equation (17)) does 
not require an underlying Keplerian disc. Hence, we can ap- 
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Figure 9. Same as in Fig. |6]but for the strong-field run 5.2-4 and after 4000 yr. The black contours enclose the regions where gas feels 
a real acceleration. Left: Centrifugal acceleration is possible only in a small part of the outfiow. Right: The general outfiow criterion fits 
the regions of gas acceleration very well. 



ply it here as well. We show the result in the left panel of 
Fig. [9] whereas the right panel shows the result of applying 
the general outflow criterion (equation (11)). For compara- 
tive purposes we also show the regions where gas feels a real, 
radial outward directed acceleration (black lines) 



dv 

ttr = -- ■ er 
at 



/ dv dxi \ _ ( ^'^ \ 
\dxi dt J ^ \dxi / 



er > 0, 



(19) 



where the time derivative is approximated by spatial deriva- 
tives. As already observed in Fig. [S] regions of acceleration 
can be found at radii > 300 AU but also closer to the z- 
axis and significantly above/below the disc. Purely centrifu- 
gal acceleration, however, only works in the outer regions 
(left panel of Fig. [9]). In the upper outflow lobe the grey 
shaded region fits the outer region of acceleration reason- 
ably well. Although in the lower outflow lobe the agreement 
is rather poor, it shows that in case of sub-Keplerian disc 
rotation the criterion clearly works better than the 30° cri- 
terion which would predict outflowing gas everywhere. Fur- 
thermore, the result suggest that the bulk of outflowing ma- 
terial, which emerges from radii > 300 AU (compare upper 
panel of Fig.jsj), is launched centrifugally also in this case. 

However, in particular in the lower outflow lobe and 
close the ^-axis, the predicted regions of purely centrifu- 
gal acceleration (grey shaded areas in the left panel) hardly 
match the regions of real acceleration (enclosed by black 
contours). Hence, as in the weak fleld case 26-4, we suppose 
that in these regions the toroidal magnetic fleld strongly 
influences the outflow dynamics. Indeed, analysing the out- 
flow with the general criterion given in equation (11) taking 
into account 5^ (right panel of Fig. [9]) shows that now the 
grey shaded regions flt the regions of real acceleration much 
better, in particular close to the z-axis. Naturally, the agree- 
ment is not perfect which is probably due to the azimuthal 
averaging process and the fact the outflow is not stationary 
as assumed for the derivation of the criterion. Nevertheless, 
it agrees remarkably well, demonstrating the importance of 
Bd). Hence, we tentatively suggest that in this run its influ- 
ence on the driving is even larger than in run 26-4 where the 
fast gas is mainly driven centrifugally. 



Next, we consider the reason for the poor collimation of 
the outflow. As mentioned before, outflow collimation is due 
to the hoop stress produced by B^. Furthermore, we again 
note that the disc has strongly sub-Keplerian rotation due 
to very efficient initial magnetic braking (see Paper I). As 
magnetic braking basically transfers a part of the rotational 
energy of the disc into energy of the toroidal magnetic fleld, 
one could naively expect in the outflow to be relatively 
strong and hence the outflow to be well collimated. How- 
ever, as already indicated in Fig. |8] is larger than Bpoi 
only in parts of the outflow suggesting a rather moderate 
toroidal fleld strength. Indeed, comparing the absolute en- 
ergy content Emag,^ stored in the toroidal magnetic fleld of 
the outflow in the strong-fleld run 5.2-4 with the weak-fleld 
run 26-4 after 4000 yr, reveals that Emag,^ is more than 3 
times larger in the latter. The reason for this is that a sig- 
niflcant fraction of the angular momentum is removed by 
magnetic braking already before the gas falls onto the disc. 
Indeed, the disc in run 5.2-4 has a 2 times lower speciflc 
angular momentum 



(20) 



Mdie 



than the disc in run 26-4. Therefore, there is less rotational 
energy left to transfer into Emag,^- Consequently, in the 
outflow and therefore also the hoop stress responsible for 
collimation are signiflcantly weaker. As at the same time 
the poloidal magnetic fleld is quite strong in run 5.2-4, colli- 
mation gets even more difficult thus resulting in the poorly 
collimated, sphere-like outflow. 



Peters et al. (2011), studying the interplay of ionising 



radiation and magnetic flelds, flnd a similar sphere-like out- 
flow structure existing over a timescale of several 10^ yr. 
The authors simulated a 1000 Mq core with a resolution of 
98 AU and a weaker initial magnetic fleld (/i = 14) than in 
our run 5.2-4. Due to the gravitational instability and sub- 
sequent fragmentation of the disc as well as the effects of the 
emerging H II regions no well deflned Keplerian disc builds 
up in their simulation. As argued above, this and the re- 
sulting weaker lead to the development of a poorly colli- 



Magnetic fields in massive star formation II 13 



t = 5000 yr 



2000 



-2000 




t = 4000 yr 



1000 



500 



-500 



-1000 



-2000 




Figure 10. Slice along the z-axis for the weak-field run 26-4 after 5000 yr (left) and the strong-field run 5.2-4 after 4000 yr (right). 
The outflow velocity Vz,out is overlaid by the contours where Vz,out equals the escape velocity Vesc (black line) and the poloidal fast 
magnetosonic velocity (white line). For the weak magnetic field case 26-4 (left), Vz,out is almost everywhere higher than -Uesc and 'Ufast- 
In contrast, for run 5.2-4 in the most parts Vz,out smaller than Vesc and exceeds Vf^st only in parts of the outflow. Note the three times 
larger spatial scale in the left panel and the different colour scaling. 



mated outflow. Hence, despite different initial conditions the 
similarities between the outflows observed by |Peters et al.| 
(2011) and us are not surprising. 



To summarise, due to magnetic braking the gas in the 
disc in run 5.2-4 is rotating relatively slowly. Therefore, 
generated by rotation is weak resulting in a poorly coUi- 
mated outflow. 



4.3 Long term evolution and jet stability 

Due to computational cost reasons we cannot follow the out- 
flow evolution over more than a few 10*^ yr. Nevertheless, we 
can try to estimate whether these outflows will persist over 
time or fall back due to the gravitational attraction of the 
central object and the ram pressure of the inf ailing gas. For 
this purpose we compare the outflow velocity Vz,out with two 
basic velocities. The first one to compare with is the escape 
velocity 



2GMsd 



(21) 



where Msd is the mass of the star + disc system. For the 
outflow to escape the gravitational potential of the central 
object Vz,out has to exceed Vesc- To take into account the 
effect of thermal and magnetic pressure in the surrounding 
gas which also might slow down the outflow, we furthermore 
compare the outflow velocity to the poloidal component of 
the fast magnetosonic velocity 



Viast 



^i,pol + Cs 



(22) 



In Fig. ^] we show the azimuthal average of outflow velocity 
Vz,out and the contours where Vz,out equals the two velocities 
mentioned above in a slice along the 2;-axis for run 26-4 (left 
panel) and run 5.2-4 (right panel). We first concentrate on 
run 26-4. As already seen in Fig.|4] below 2000 AU the high- 



est velocities occur close to the z-axis. Furthermore, it can 
be seen that the outflow velocity exceeds the escape velocity 
and also the fast magnetosonic velocity in most parts of the 
outflow. 

In Fig. [TT] we plot the ratio of Vz,out to 
■^^= and va (equation (js )) for run 26-4 along a vertical line 
at a radius of 100 AU ana 300 AU. As before the quantities 
are averaged azimuthally. The outflow velocity exceeds the 
escape speed by up to one order of magnitude (bottom panel 
of Fig. 11 ) and the poloidal Alfvenic velocity even by up to 
two orders of magnitude although the average ratio is 5 - 



10 (middle panel of Fig. 11). We attribute the higher value 
of ^^z,out/'^^A,poi at a radius of 300 AU to somewhat lower 
densities in this region (compare bottom panel of Fig. |4|. 
The absolute outflow velocity, however, decreases with in- 
creasing radius as already seen in Fig. [4] and [To] Interest- 
ingly, over the whole outflow extension the gas velocity is 
comparable to or only slightly larger than va- This possibly 
points to some kind of self-regulation in the outflow where 
the gas speed is held in the trans- Alfvenic range by inter- 
nal shocks (see Fig. [2]). Once the gas reaches a significantly 
super- Alfvenic speed the flow cannot be stabilised any more 
leading to instabilities and shocks which in turn reduce the 
velocity to trans- Alfvenic speeds (see Ouy ed et al.|2003 , and 
upper panel of Fig. 11 ). Beyond the shock the gas gets reac- 



celerated by magneto-centrifugal forces (Section 4.1.2) until 
it shocks again. Interestingly, the gas speed is comparable 
to Va and not to VA,poi, i.e. has to be taken into ac- 
count. This indicates that Bpo\ is not the main agent stabil- 
ising the flow although it certainly contributes ( Ray||198T ). 
Indeed, Appl & Camenzind (1992) show that jets with a 



toroidal magnetic field are even more stable than jets with a 
purely poloidal field. In our case we therefore suppose that 
the stability of the jet is significantly enhanced by B^. 

Despite ongoing fragmentation of the protostellar disc 



(see Paper I for details) , the driving of the outflow seems not 
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Figure 11. Ratio of Vz,out to the Alfvenic speed va (top), to the 
poloidal Alfvenic speed ^'A,pol (middle) and to the escape speed 
'^esc (bottom) for run 26-4 along a vertical line at a radius of 100 
AU (black hue) and 300 AU (red line). The outflow velocity is 
signiflcantly higher than i^esc 

and 'UA,pol whereas it is comparable 
to Va, probably indicating a self-regulated outflow speed. 

to decline over time. As furthermore Vz,out is significantly 
larger than Vesc and Vfast , it can be expected that the outflow 
in run 26-4 will escape the gravitational potential of the 
central star/disc-system and finally will leave the core even 
when taking into account the ram pressure of the infalling 
gas. 

Considering run 5.2-4 in the right panel of Fig. [lO] re- 
veals a markedly different situation. Comparing Vz,out to 
Viast shows that the gas is moving super- Alfvenic only in 
parts of the outer wings of the outflow. Furthermore, a com- 
parison with ^;esc shows that the outflowing gas is moving 
almost everywhere with velocities below the escape veloc- 
ity. In general, Vz,out exceeds Vesc and Vfast by a factor of at 
most 2. Therefore, when taking into account the additional 
deceleration of the outflow by the ram pressure of the in- 
falling material, it might be possible that the outflow will 
re-collapse to the star/disc-system. 

To summarise, the outflows of the two simulations pre- 
sented so far do not only differ in their morphology and their 
kinematics but possibly also in their longer term evolution. 



t = 4000 yr v, [km s'] 
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Figure 12. Slice along the z-axis in run 10-20 after 4000 yr show- 
ing the outflow velocity. Two distinct outflow components are 
clearly visible, the inner, well collimated, fast jet and the outer, 
low-velocity outflow. 

The well collimated outflow is likely to overcome the gravi- 
tational potential of the central protostar and leave the core 
whereas the sphere-like, slowly expanding outflow might be 
only a short-lived, transient feature in the very early phase 
of massive star formation possibly recoUapsing again. 

4.4 The influence of the initial conditions 

So far only the outflows in the runs 26-4 and 5.2-4 have 
been considered in detail. We have chosen these particular 
runs as they give representative examples for the outflows 
observed in the other simulations. The remaining outflows 
reveal qualitative similarities to one of the two outflows pre- 
sented before. One of the main results of the previous sec- 
tions is that the different morphologies can be attributed to 
the different velocity structure in the disc resulting in dif- 
ferent radial positions where the outflows are launched and 
variations in the strength of the hoop stress (B^) responsible 
for outflow collimation. To further confirm this, in the fol- 
lowing we consider the whole set of simulations and connect 
the outflow properties to the properties of the protostellar 
disc. 

We start with the remaining weakly magnetised runs. 
The outflows in the runs 26-20, 26-0.4 and 26-0.04 show 
similar morphologies as the outflow in run 26-4. They are 
reasonably well collimated with collimation factors between 
2 in run 26-20 and ^ 4.5 in run 26-0.4 and have outflow ve- 
locities preferentially parallel to the 2;- axis. Interestingly, in 
the runs 10-20 and 5.2-20, whose outflows are rather poorly 
collimated, a fast and well collimated jet close to the z-axis 
develops after about 2500 yr embedded in the slowly ex- 
panding, poorly collimated outflow. Similar applies to run 
26-20 with the difference that by the end of the simulation 
the initial, slowly expanding, poorly collimated outflow is 
overtaken by the faster, more collimated jet. For demonstra- 
tive purposes, in Fig. [12] we show the situation in run 10-20 
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after 4000 yr. We mention that similar outflows consisting 
of two distinct components, an outer, slowly expanding and 
an inner, fast component have recently also been observed 
in low-mass star formatio n simulations (Banerjee Pudritz 



2006] IMachida et al. 120081 |Hennebelle Fromang|2008| |To- 
mida et al.|2010| [Puffin fc Pudritz|2009| [Puffin et al.|2C)ll| ). 



Analysing the magnetic field line structure in the runs 
mentioned before reveals situations similar to that in run 
26-4. In the inner region close to the symmetry axis of 
the jet, where the gas gets accelerated, the magnetic field 
is only weakly dominated by the toroidal field component, 
i.e. B(f)/Bpoi ^ 10. Furthermore, the field lines at the disc 
surface are inclined by more than 30° w.r.t. to the verti- 
cal axis. As the discs rotate with Keplerian velocities just 
like in run 26-4 or as in run 5.2-20 with at least about half 
the Keplerian speed, we therefore suppose that the gas gets 
launched from the discs by centrifugal acceleration. Apply- 



ing our outfiow criterion (equation ( 17 )) for these runs shows 
that above the disc the region where purely centrifugal ac- 
celeration works is not as extended as in run 26-4 (Fig. |6|. 
With the general criterion for magneto-centrifugal accelera- 
tion (equation ( 11 )), however, we can fit the outfiow regions 
by far better. Hence, considering the jets as purely centrifu- 
gally driven might not be an appropriate choice too far from 
the disc. The evidence shows that the pressure gradient of 
the toroidal magnetic field contributes significantly to the 
gas acceleration although, as mentioned before, the launch- 
ing from the disc itself is most likely to due to centrifugal 
acceleration. The application of the two criteria also suggest 
that the dynamics of the outer, slowly expanding outfiows 
in the runs 26-20, 10-20 and 5.2-20 are mainly determined 
by B^. 

All remaining runs show poorly coUimated, low veloc- 
ity outfiows with outfiowing gas emerging almost radially at 
radii > 100 AU and gas infall close to the z-axis. Further- 
more, the magnetic field line properties in the centre are 
similar to run 5.2-4 and the protostellar discs driving the 
outfiows are all clearly sub-Keplerian. Moreover, applying 
the outfiow criteria derived in Section [3] gives very similar 
results as for run 5.2. Therefore we conclude that in these 
runs the outfiow driving mechanism is very similar to that 
in run 5.2-4 , i.e. gas is launched centrifugally from the disc 
at large radii (r > 100 AU) with support by i?^. At larger 
distances is mainly responsible for the further acceler- 
ation (see Section 4.2.2[ ). To remind the reader, the poor 
collimation is a consequence of the weak toroidal magnetic 
field in the outfiow and of the large launching radii. The 
larger launching radii are a result of the fact that close to 
the centre gravity is too strong to be overcome by the weak 
centrifugal force. 

In Fig. [13] we show the dependence of the outfiow mor- 
phology on the initial conditions. It can be seen that only 
for the runs with weak magnetic fields (/x ^ 26) or high 
rotational energies (/3rot = 0.20) do well collimated jets oc- 
cur. Moreover, all runs with strongly sub-Keplerian discs 
have poorly collimated outfiows (see Fig. 8 in Paper I|). This 



clearly demonstrates that for a collimated, high velocity jet 
to form a disc with (almost) Keplerian rotation is necessary. 
In contrast, if the disc is rotating strongly sub-Keplerian, 
both outfiow speed and collimation decrease significantly. 

Next, we consider the infiuence of the initial conditions 
on global outfiow properties like mass and momentum. The 
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Figure 13. Phase diagram of magnetic field and rotational en- 
ergy (see Table Q showing the results of the simulations concern- 
ing the question of jet formation. For runs with weak magnetic 
field (/i = 26) or high rotational energies well collimated, cen- 
trifugally driven jets develop (red squares). In contrast, for the 
remaining runs only slowly expanding, poorly collimated outflows 
are generated (black circles). 



results for each simulation are listed in Table (2] For the cal- 
culation of mass and momentum we only take into account 
gas more than 47 AU above or below the midplane. By doing 
so, we avoid the outfiow quantities to be affected by mate- 
rial in the protostellar disc. For the momentum we only con- 
sider the z-component as this is its preferred direction. The 
outfiow masses are of the order of 0.1 - 1 M(;) s ignificantly 
lower than the sink masses (see also Section 4.5). With mo- 
menta of the order of 0.1 - 1 Mq km s~^ we obtain average 
velocities of a few km s~^ lower than typically found in ob- 
servations (e.g. [Beuther et al.|2002b| see also Section 5.3 for 
a detailed discussion). 

In general, for a fixed magnetic field strength both cen- 
trifugal and Bd) dominated acceleration predict that the 
power of the outfiow should increase with the rotation speed 
of the protostellar disc. Thus, for larger amounts of initial 
rotational energy and hence faster disc rotation, the outfiow 
mass and momentum are expected be higher as well. Indeed, 
comparing runs with equal /j, in Table |2] reveals that except 
for run 26-20 the outfiows are in general more powerful for 
higher /3rot- We suppose that in run 26-20 the generation of 
the outfiow might be slightly delayed as the material falling 
onto the disc has a high excess of angular momentum which 
has to be removed before moving further inwards. Therefore, 
the accretion rate and consequently also the outfiow rates 
Mout and Pout are lower. 

Interestingly, there is no clear trend recognizable in 
Mout and Pout when considering simulations with fixed Prot 
but varying /i. In general, for runs with /i ^ 5.2 there seem 
to be only small variations within a factor of about 5 in Mout 
and Pout (see TablejS]). We note that we exclude run 26-20 in 
this consideration as it does not fit in this trend. For strong 
magnetic fields (/i = 2.6), however, there seems to be a rapid 
decline in outfiow power. Naively, it could be expected that 
the power of the outfiows would increase with magnetic field 
strength. This is not the case since for stronger magnetic 
fields disc rotation gets more and more sub-Keplerian (see 
[Paper I[ for a detailed discussion). This results in a weaker 
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centrifugal acceleration and a weaker toroidal magnetic field 
responsible for further accelerating the outflow. Therefore, 
for both driving mechanisms the outflow power is expected 
to decrease for strongly sub-Keplerian discs in agreement 
with our observations. However, we emphasise that, as the 
outflow morphologies differ significantly between the indi- 
vidual simulations, the results of this comparison have to be 
taken with care. 

In summary, we find that the generation of a fast 
(vz,out — 10 km s~^), well collimated jet depends on 
the build-up of a Keplerian disc. In these jets the gas is 
most likely launched centrifugally from the disc. Somewhat 
above/below the disc, however, we expect the Lorentz force 
to come into play contributing to the acceleration of the 
gas. In contrast, for sub-Keplerian discs slowly expanding, 
sphere-like outflows develop in which the gas is launched 
centrifugally from the disc as well although at larger radii. 
Also here 5^ quickly plays an important role in accelerating 
the outflow further. Although the morphologies are very dif- 
ferent, the global properties like mass and momentum vary 
only within a factor of 5 showing a decrease towards high 
magnetic field strengths. 



4.5 The ejection/accretion ratio 

Next, we compare the mass outflow rates Mout with the mass 
infall rates Min- To begin with, we limit our consideration 
to the fiducial runs 26-4 and 5.2-4. In Fig. ^] we show the 
ratio of outflow to infall rate Mout / Min as a function of the 
distance d from the centre for different times. The outflow 
and infall rates are computed on the surface of a cylinder of 
radius d and height 2d centred around the midpoint of the 
simulation domairj^ In run 26-4 Mout /Min is not constant 
over the outflow extension but varies between 0.1 and 0.4 for 
d > 500 AU. We attribute these variations to the internal 
shocks occurring in the outflow (compare Fig. |2] and 10). 
Although being somewhat higher, in general the observed 
range roughly agrees with theoretical estimates of ^ 0.1 ob- 



tained for the centrifugal acceleration mechanism ( Pudritz 
Norman||l986| [Pelletier Pudri"t^[r992J . The high value 
of Mout / Min at d < 500 AU is due to the circulation flow in 
the outer parts of the outflow (r ~ 200 - 300 AU, see bot- 
tom panel of Fig. |5|. As the outflow speed of this material 
is not sufficient to escape the gravitational potential, it falls 
back to the disc thus not contributing to the outflow rate 
at larger distances. The shape of Mout/Mm in run 5.2-4 is 
qualitatively very different from that in run 26-4. The ratio 
seems to saturate around a value of ~ 0.1 for late times and 
small distances (d ^ 400 AU) in agreement with theoretical 
estimates. For larger distances, however, Mout /Min increases 
and reaches values up to 0.6. However, also for this run the 
measured Mout /Min is in reasonable agreement with theo- 
retical predictions for centrifugal acceleration even though 
partly somewhat too high. 

Next, we consider Mout /Min for the whole set of simula- 
tions performed. To avoid the problem that even for an indi- 
vidual outflow this ratio is not constant in time and position, 
we follow a simpler approach than before. For this purpose 



^ We have chosen a cylinder instead of a sphere to lower numerical 
errors when computing the flux through the surface. 



we approximate Mout /Min by calculating the ratio of the 
total outflow mass to the total mass accreted onto the sink 
particles (see Table [2]). The result is shown in Fig. 15 Ex- 
cluding one exemption (run 10-20) the ratios are all smaller 
than 0.5. The values have a mean of ~ 0.3 which agrees 
very well with other numerical work ( Tomisaka||1998 2002 



Hennebelle Froma ng 2008 : [ Puffin et al.|2011t an d is only 
slightly higher than analytical estimates of ~ 0.1 (Pudritz 
&: Norman 1986) [Pelletier k Pudritz]|1992t . However, the 
scatter of the individual values is significant ranging from 
about 0.02 in run 26-0.04 up to 0.9 in run 10-20. Interest- 
ingly, run 26-20 has a value of Mout/Min around 0.3 in ac- 
cordance with the other runs although the absolute values 
in this simulation are very small (see Table|2]). This conffims 
the explanation given in Section |4.4| that in run 26-20 the 
whole outflow ejection process is slowed down due to the 
high excess of angular momentum which causes lower accre- 
tion rates and in turn also lower outflow rates. Moreover, 
Mout/Min is generally higher for higher /3rot- As discussed 
in Section |4.4| this is a consequence of the fact that the 
outflow power increases with firot whereas at the same time 
the accretion rates decrease due to an enhanced centrifu- 
gal support. For varying /x, however, there is no clear trend 
in Mout / Min recognizable except a possible drop-off for the 
highly magnetised runs with /i = 2.6 which was already ob- 
served in Section [X4l for the absolute outflow quantities. 

Although the derived values of Mout / Min are in general 
somewhat higher than analytical predictions, they agree rea- 
sonably well with observational results ranging from about 
0.1 to 0.3 ([Richer et al.|2000||Beuther et al.|2002"bl|Klaassen 



et al. '2011). Furthermore, a ratio of ^ 0.3 is often used 



to estimate accretion rates of massive protostar from ob- 
served outflow mass rates (Beuther et al. 2002a|b 2003 



2004 



Zhang et al.[[2005 ). As this value fits reasonably well 



over a wide range of initial conditions, we suppose that ac- 
cretion rates derived that way seem to be reliable. 



5 DISCUSSION 



5.1 Launching mechanisms 



Our analysis presented in the Sections [4.1.2[ and [4.2.2[ sug- 
gests that in all simulations the gas is launched from the 
disc by centrifugal acceleration (BP82). At larger distances 
above/below the disc the pressure of the toroidal magnetic 
field significantly contributes to the acceleration of the gas. 
Although Lynd en-Bell[ ([1996| [2003| describe outflows whose 
dynamics are determined by the toroidal magnetic field, we 
emphasise that our outflows are for the most parts (in time 
and space) not such so-called magnetic tower flows. In fact, 
the effect of is also implicitly contained in the MHD wind 



theory (equation (11)). Therefore, it is more appropriate to 
consider our outflows in the framework of MHD wind theory 
as it contains both centrifugal acceleration and the magnetic 
pressure (B^) effects. In such a magneto-centrifugally driven 
wind the gas is ejected from the disc surface and then gets 
accelerated towards larger heights. In contrast, in a mag- 



netic tower (Lynden-Bell 1996 2003.) gas gets accelerated 



only at the tip of the outflow - there is no sustained mass 
loss coming from the disc. Tower flows are transient phases 
that arise only in the earliest stages of outflows. 
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Figure 14. Ratio of the mass outflow rate over infall rate in the weak-fleld run 26-4 (left) and the strong-fleld run 5.2-4 (right) as a 
function of the distance from the centre. The values in both runs agree with theoretical estimates and other numerical work. 
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Figure 15. Ratio of outflow mass to sink mass at the end of each 
simulation as a function of /j, (see Table Q . Equal symbols de- 
note equal initial rotational energies. The calculated values scat- 
ter around a mean of ~ 0.3 in rough agreement with theoretical 
estimates. 



For simulations with weak initial magnetic fields or high 
rotational energies the situation is as follows: shortly after 
the formation of the first sink particle, a sub-Keplerian disc 
builds up driving a slow (^ 1 km s~^) outflow whose driving 
is probably dominated by the magnetic pressure. As soon 
as a well defined large-scale Keplerian disc develops, in the 
inner region a fast (~ 10 km s~^), well collimated jet is 
launched by centrifugal acceleration. At larger heights and 
radii > 200 AU, however, significantly contributes to the 
acceleration. An outer, slowly expanding outflow component 
is still clearly visible at the end of the simulations 26-20, 



10-20 and 5.2-20 (see Fig. 12). Such outflows consisting of 
two components, an outer, slowly expanding and an inner, 
fast component have also been observed in simulations of 
low-mass star formation ( Banerjee Pudritz|2006||Machida 



et al.|2008||H ennebelle Fromang 2008! |Tomida et al.|2010[ 

J. 



Duflfin & Pudritz 2009; Duffin et al. 201 



For runs with strong magnetic fields, the protostellar 
discs are sub-Keplerian all the time due to strong magnetic 
braking (see |Paper I| for a detailed discussion). Therefore, 
during the time covered by our simulations no fast, central 



jet component develops and only a slowly expanding, poorly 
collimated bubble is present. The gas seems to be mainly 
launched centrifugally although the toroidal magnetic field 
has probably a large impact on driving the outflow already 
close to the disc. We tried to clarify the situation by ap- 
plying a criterion derived in Section [3] depending on the 
actual position, the protostellar mass and the local rotation 



speed (equation (17)). In general, using the criterion for the 
present simulations indicates that it is applicable despite the 
constraints like axis-symmetry, stationarity and corotation 
assumed for its derivation. We show that for strongly sub- 
Keplerian discs it acceptably fits the region where the out- 
flow gets launched (see left panel of Fig.[9| in contrast to the 
[Blandford &: Payne] criterion which fails when analysing the 
sub-Keplerian cases. Together with the results from run 26-4 
this demonstrates that simply considering the inclination of 
the magnetic field w.r.t. the z-axis might not be sufficient to 
determine the driving mechanism nor the region where the 
outflow gets launched. On the one hand even for inclination 
angles smaller than 30° the outflow can be driven centrifu- 
gally (see Fig. [s] and |6| whereas on the other hand there 
are situations where purely centrifugal acceleration does not 
work despite inclination angles larger than 30° (see Fig. [s] 
and [9]). 

Again, we emphasise that in all runs the mechanism 
responsible for outflow launching is centrifugal acceleration. 
Only at larger heights do the effects of the toroidal magnetic 
field come into play. The situation is therefore not such that 
the outflow is either driven centrifugally or by the toroidal 
magnetic pressure. Such a classification often found in lit- 
erature is an oversimplification not describing the situation 
appropriately as both processes play an important role. We 
have shown that there is a continuous transition, from cen- 
trifugal dominated acceleration to acceleration dominated 
by B(f). These two regimes are both implied in the solution 
of the stationary MHD equations presented in Section [3] (see 
also Section 4.1 of Spruit || 1996 ). The different outflow mor- 
phologies mainly result from the varying strength of re- 
sponsible for collimation as for strongly sub-Keplerian cases 
B(f) is too weak to properly collimate the outflows. 

As mentioned before, several authors have simulated 
outflows consisting of two different components. However, 
their conclusions as which component is driven by which 
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mechanism partly contradict each other. On the one hand, 
several authors ( Banerjee &: Pudritz|2006 Duffin & Pudritz 



2009) claim that the outer, larger scale outflow is driven 
by the magnetic pressure (Lynden-Bell 2003 ) whereas the 
inner, fast outflow embedded in the first one is considered 



as a centrifugally driven jet (BP82). On the other hand. 



Machida et al. (2008) and Tomida et al. (2010) find the sit- 



uation to be exactly vice versa. All authors argue on the 
basis of an analysis of the field line structure such as the 
ratio of to Bpoi or the field line inclination w.r.t. the 
protostellar disc. As we argued above, this kind of consider- 
ation alone can be misleading which might have caused the 
contradictory results. Hence, it would be interesting to apply 
the criterion derived in this work to confirm the results in 
an independent way. This was recently done by [Duffin et al.| 
(2011) in simulations of low-mass cores. The results sup- 



port our interpretation of a magnetic pressure driven, outer 
outflow component and a centrifugally driven, inner com- 
ponent. Furthermore, this demonstrates the applicability of 
our criterion in equal measure to low-mass and high-mass 
protostellar outflows without any restrictions. 

Hennebelle Fromang] ([2008) studying the influence 



of the initial magnetic field strength on the formation of 
outfiows around low-mass stars find striking similarities to 
our results, i.e. a well collimated, fast outfiow component 
embedded in a slowly expanding component for weak mag- 
netic fields and a poorly collimated outfiow for strong fields. 
Their interpretation of a centrifugally driven wind for the 
latter case and a magnetic tower fiow for the former case, 
however, is in opposition to our interpretation. Hence, it 
would be interesting to apply our criterion derived here to 
confirm this conclusion in another way. Differences might 
also result from the fact that Hennebelle & Fromang ( 2008 ) 



consider a low-mass core which is less gravitationally unsta- 
ble. Therefore, in particular in the strongly magnetised case 
the situation might be less sub-Keplerian possibly leading 
to different results. 

In general, our results agree remarkable well with those 
of low-mass star formation simulations even over the wide 
range of initial conditions covered in out work. However, the 
good agreement between our and the low-mass case results 
could not be expected a priori as in contrast to low-mass 
cores the cores presented here are highly gravitationally un- 
stable (about 56 Jeans masses). Therefore the dynamics and 
timescales of the whole system can differ significantly from 
low-mass cores. Another difference, though not covered in 
this work, is the infiuence of initial turbulence and radia- 
tive feedback of the protostars. However, as we are still in a 
very early evolutionary stage, we suppose that the effect of 
radiation is still rather limited. In contrast, the infiuence of 
turbulence which we will examine in a subsequent paper is 
expected to be significant (see also [Hennebelle et al.||2011| ) . 



5.2 The evolution and impact of outflows 

For typical massive star forming regions the observed mass- 
to-fiux ratios are usually only slightly supercritical, /x < 



5 (Falgarone et al.p008i Girart et al.||2009 Beuther et al. 



2010). Based on the results obtained in this work, the ques- 



tion arises how well collimated, fast outfiows frequently ob- 
served around high-mass protostellar objects can form (e.g. 
Beuther et al.||2002bl 120041 iBeltran et al.||2011|). As we 



showed that the formation of such outfiows is linked to the 
build-up of Keplerian discs, the question reduces to how the 
magnetic braking efficiency can be reduced allowing Keple- 
rian discs to form. A possible solution is provided by the 
poorly collimated outfiows themselves evacuating the region 
above and below the disc therefore significantly reducing 



the magnetic braking efficiency ( [Mouschovias Paleologou 
Furthermore 



1980) 



misalignment between the rotation 



axis and the initial magnetic field ( .Hennebelle Ciardi 



2009 



Ciardi k Hennebene|[20Tot or non-ideal MHD effects 



like ambipolar diffusion, ohmic dissipation or the Hall effect 
could reduce the magnetic braking efficiency, in particular 
in later stages (e.g. |Dapp k Basu 2010| |Li et al.||201l| . Ah 
these possibilities lead to a successive build-up of Keplerian 
discs during the protostellar evolution towards Class I/II 
objects, hence allowing the generation of well collimated, 
high- velocity jets (see also PaperT{ for a more detailed dis- 
cussion) . 

As shown in Fig. [s] the position- velocity (PV) diagram 
shows a typical "Hubble law "-like behaviour also seen in 
observations (e.g. |Lada Fich| |1996^ 'Beltran et al.| |2011 ). 
Such an behaviour is a natural outcome of several wind mod- 
els like the jet-driven bow shock model or the wind-driven 



shell model (see Cabrit et al.||1997^ for an overview). How- 



ever, the comparison between our results and those models 
is not straightforward as in those wind models the veloc- 
ity increase with distance is only an apparent and no real 
acceleration. This apparent acceleration is a consequence of 
the inclination of the outfiow w.r.t. the plane of the sky and 
the fact that mainly gas swept up in the bow shock and not 
inside the outfiow lobe contributes to the observable emis- 
sion. In contrast, in our case the observed acceleration is real 
and is mainly due to the gas inside the bow shock structure 
whereas the material in the shock itself has mainly negative 
velocities (compare Fig. [4]). Hence, an observed acceleration 
of the bulk velocity in outfiows can also be real and not only 
due to projection effect. 

However, an interesting point made by those wind mod- 
els is the fact that the outfiows partly consist of gas en- 
trained from the ambient medium. We tested this statement 
by implementing a passive mass scalar in run 26-4 which re- 
veals an outfiow morphology similar to that assumed in the 
wind models. This mass scalar is set to 1 in a disc of 75 
AU above and below the midplane and gets advected with 
the gas. Hence, we can test how much of the gas in the out- 
fiow stems from this disc and how much from the ambient 
medium. The analysis shows that by the end of the simula- 
tion about 92% of the gas mass in the outfiow stems from 
the disc. Hence, most of the gas gets accelerated from the 
bottom of the outfiow up to the tip and only a small part 
is entrained from the ambient medium. A possible expla- 
nation might be the high infall velocities which are a result 
of the highly gravitationally unstable configuration. We sug- 
gest that the gas hitting the bow shock it is mainly defiected 
sidewards and gets channelled downwards along the shock 
rather than getting stuck and being entrained with the out- 
fiow. The situation might change if the core is less gravi- 
tationally unstable and the infall velocities are smaller. We 
tentatively suppose that in this case it is more likely that 
the gas gets entrained by the outfiow. 



Yorke k Sonnhalter (2002) and Krumholz et al. (2005) 



suppose that outfiows create low density regions along the 



Magnetic fields in massive star formation II 19 



polar direction which aid the formation of massive stars by 
channelhng the radiation outwards, therefore diminishing 
the radiation pressure and allowing accretion to proceed. 
As seen in bottom panel of Fig. [5] the material ejected 
from the disc indeed carves a low density region into the 
vicinity. It would be interesting to see whether the density 
decrease in our simulations is sufficient to allow radiation 
produced by the protostar to escape efficiently as recently 
seen in radiation-hydrodynamical simulations by Cunning-] 
ham et al. (2011 ). Furthermore it would be of interest to see 



to what extent the radiation pressure would decollimate the 
outflow as proposed recently by |Vaidya et al. (2011) on the 
basis of 2-dimensional jet simulations. 



5.3 Comparison to observations 

In Tabled we have listed the mass and momentum outflow 
rates of all runs performed. The calculated mass-loss rates 
vary between 1.6 • 10"^ and 2.7 • 10"^ Mq yr"^ This is in 
agreement with typically observed mass loss rates in high- 
mass star forming regions ranging from 10~ ^ up to a few 



10 ^ M(;) yr ^ although at the lower end (e.g. Beuther et al. 



2002b| [Zhang et al.| 2005] [Wang et al.|2011||Ren et al.|2011| ). 



Like our simulations those observations refer to the very 
early stage of massive star formation although the calculated 
dynamical timescales are of the order of a few 10^ up to a 
few lO'^ yr, typically somewhat larger than ours. However, 
as for a large fraction of the observed systems no signs of 
H ii-regions have been detected, the outflows are probably 
still magneto-centrifugally driven winds which is why our 
mass- loss rates are comparable. 

The calculated momentum outflow rates in our models 
range from 4.0 • 10~^ to 4.7 • lO"'^ Mq km s"^ yr"\ This 
is at the lower end usually reported by observations of mas- 



to 10~^ M(7) km s~^ yr"^ (e.g. Beuther et al. 


2002b 


Zhang 


et al.||2005t Shi et al.||2010 


Wang et al.||2011 


). As the out- 



flow velocity in centrifugally driven jets is coupled to the 
rotation velocity (Michel "1969J [Pelletier Sz Pudritz'> 1992| ), 
the maximum outflow velocities which can be reached in our 
simulations, are of the order of 10 km s~^ (see equation ([2])), 
about one order of magnitude lower than typically observed. 
Therefore, the lower values of Pout in our runs are not sur- 
prising. However, we expect our momentum outflow rates to 
increase and to become comparable to observations if we run 
the simulations over comparable timescales of a few 10^ yr 
and with a higher resolution which allows for faster outflow 
material arising from deeper in the gravitational potential 
well. 

Observational results on outflow morphologies are still 
somewhat ambiguous. One the one hand, outflows around 
massive protostars are usually found to be less collimated 
than their low mass counterparts revealing coUimation fac- 
tors of the order 



/u| 
f 



of 1 - 2 (e.g [Ridge Moore 2001; Wu 
et al.|2004t . On the other hand, [Beuther et al.[ ( |2002a,,200^ 
also find highly collimated, massive outflows with coUima- 
tion factors as high as 10. Based on this, 'B euther Shep-] 
herd (2005) propose an evolutionary scenario in which well 



collimated, magnetically driven outflows occur in the very 
early phase of massive star formation. In their further evo- 
lution the outflows get progressively less collimated due to 
the build-up of H ii-regions. According to this, early-stage 



outflows as presented here should all be more or less well col- 
limated. Obviously this is not the case, as only runs with a 
weak magnetic field (/i = 26) reveal well collimated outflows 
with collimation factors up to 4.5. Therefore, our results sug- 
gest that the collimation not only depends on the evolution- 
ary stage but also on the initial conditions of the molecular 
cloud core, in particular on the magnetic field strength. 



5.4 The evolution of outflow collimation 

Our results suggest that during the very early stages (10^ 
- 10^ yr) outflows in typically magnetised, massive cores 
(/i < 5) should rather be poorly collimated with collimation 
factors of 1 - 2 instead of 5 - 10. This agrees with a number of 
observations of outflows around young massive protostellar 
objects (e.g. [Ridge k Moore|200ll [Torrelles et al.|2003|[Wu 



et al.|2004[ [Sollins et al.|2004|[Surcis et al.|2011t . Therefore, 
we suggest that in the earliest stage, i.e. even before the 



scenario described by Beuther & Shepherd (2005) applies. 



outflows are indeed rather poorly collimated except in case 
of an unusually weak magnetic field. In their further evo- 
lution, however, the collimation will increase quickly due to 
the development of a fast, central jet coupled to the build-up 
of a Keplerian disc. Such a behaviour can in particular be 
seen in the runs 26-20, 10-20 and 5.2-20. Hence we suggest 
that, beside tracing later stages as proposed by [Beuther fc| 
Shepherd (2005), poorly collimated outflows around massive 
protostars could also trace the very early stage of massive 
star formation. Hence, one might have to be careful when 
trying to infer the stage of star formation from the collima- 
tion of the outflow. 

To confirm such an evolutionary scenario as proposed 
here, observations of the earliest stage of massive star forma- 
tion would be necessary in combination with magnetic field 
measurements. Unfortunately, to date such observations are 
difficult to obtain and therefore rather rare. However, there 
is an interesting observation which supports the picture of 
very early stage, poorly collimated outflows successively col- 
limating over time. Observing two spatially adjacent, mas- 
sive protostars in the star forming region W75N, Torrelles"] 



et al.|(2003) and Surcis et al. (2011) find the younger of the 



two having a spherical outflow whereas the more evolved 
protostar has a well collimated outflow. Due to the close 
proximity to each other, the authors expect the environmen- 
tal conditions to be very similar and hence not to cause the 
different morphologies. Therefore the authors conclude that 
the differences are rather a consequence of different evolu- 
tionary stages and that the younger, poorly collimated out- 
flow is possibly only a transient feature. This interpretation 
fits perfectly into our evolutionary scenario in which after 
an initial phase a poorly collimated outflow gets overtaken 
by a well collimated jet. 



6 CONCLUSIONS 

We have studied the collapse of massive molecular cloud 
cores with varying initial rotational and magnetic energies. 
The mass-to-flux ratios of the magnetically supercritical 
cores range from 2.6 up to 26. The cores have rotational en- 
ergies well below the gravitational energy and contain about 
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56 Jeans masses. Hence, they are highly gravitationally un- 
stable and possible sites of massive star formation. In this 
work we focussed on the launching mechanism and the prop- 
erties of early outflows and jets. Furthermore, based on the 
stationary, axisymmetric MHD equations we derived a gen- 
eralised criterion to determine the launching mechanism of 
the outflows. The criterion is applicable to the whole outflow 
region and to situations with sub-Keplerian disc rotation. In 
the following we summarise our main findings. 
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and 



1. We show that our outfiow criterion (equations 
17) can successfully distinguish between the different driving 
mechanisms and that it works over the complete extension 
of the outfiows. We discuss an example where the outfiow 
is centrifugally driven up to about 800 AU thereby demon- 
strating that considering only the ratio of to Bpo\ is not 
sufiicient to determine the driving mechanism. Furthermore, 
we successfully apply the criterion to runs with strongly sub- 
Keplerian protostellar discs where the frequently used 30°- 
criterion for the inclination of magnetic field lines fails. 

2. Depending on the initial magnetic field strength two 
types of outfiows are observed: Well coUimated, high- 
velocity jets for runs with weak magnetic fields (/i ^ 26) and 
high rotational energies (Prot — 0.20) and poorly collimated, 
slowly expanding outfiows for runs with strong fields. The 
development of fast jets is coupled to the build-up of Kep- 
lerian discs. In none of the strongly magnetised runs with 
strongly sub-Keplerian disc rotation a well collimated jet is 
observed. 

3. In all runs centrifugal acceleration is responsible for 
launching the gas from the discs. With increasing distance 
from the disc the pressure gradient of the toroidal magnetic 
field progressively contributes to the further acceleration of 
the gas. For runs with sub-Keplerian discs the outfiows are 
centrifugally dominated at large radii (> 100 AU) where 
gravity is sufficiently reduced to be overcome by the cen- 
trifugal force. In these outfiows seems to play a more 
important role in the driving than in the fast jets. 

4. The morphological differences of the outfiows are 
mainly due to the varying strength of the hoop stress re- 
sponsible for outfiow coUimation. For sub-Keplerian discs 
the generation of happens more slowly, so that their cor- 
responding outfiows are less collimated. 

5. We show that an outfiow can be maintained despite the 
fragmentation of the protostellar disc. Furthermore, knotty 
outfiow structures can also be produced in continuously fed 
jets by gas repeatedly experiencing shocks and reacceler- 
ating again even far from the disc. For such outfiows the 
toroidal magnetic field seems to contribute significantly to 
the overall stabilisation. 

6. The observed mass and momentum outfiow rates are 
of the order of 10"^ Mq yr"^ and lO"'^ Mq km s"^ yr"\ 
respectively, in reasonable agreement with observational re- 
sults. The calculated mass ejection/accretion ratios scatter 
around a mean of 0.3 in agreement with both theoretical 
estimates and observational results. 

7. Based on the results of the strongly magnetised simu- 
lations (/i ^ 5), we suggest an evolutionary scenario where a 
poorly collimated outfiow is typical for the very early stage 
of massive star formation. Over time the outfiow coUimation 
will increase due to the development of a well collimated, fast 
jet overtaking the slowly expanding outfiow. This picture is 



also supported by observations. Furthermore, analysing the 
sphere-like, slowly expanding outfiows suggest that they are 
possibly only transient features which might re-collapse dur- 
ing their further evolution. 

We showed that a simple approach such as the value of 
B(j) I Bpo\ or the inclination of magnetic field lines is not suf- 
ficient to determine the launching mechanism of outfiows. 
This could be the reason for the partly confiicting results 
found in literature. We therefore strongly suggest to use a 
self-consistent criterion which is applicable to the whole out- 
fiow region as well as sub-Keplerian discs frequently found in 
numerical simulations (e.g. Mellon Li|2008] Hennebelle & 



Fromang|20"Q8l [Dufiin k Pudritz|2009| [Seifried et al.|2011| ). 

Moreover, we emphasise that a separation in purely cen- 
trifugally or magnetic pressure driven winds is probably an 
oversimplification not describing the situation appropriately. 
We showed that magneto- centrifugal acceleration has two 
regimes, a centrifugally dominated and a 5^ dominated one 
and that within a realistic outfiow there is a continuous tran- 
sition from one regime to the other. 

A growing number of observations of discs and well col- 
limated, bipolar outfiows around high-mass protostars (see 
[Beuther Shepherd| |2005l [Cesaroni et~ar] |2507l for recent 
reviews) support a high-mass star formation scenario via 
disc accretion similar to low-mass star formation. On the 
other hand, observations also show that prestellar cores with 
masses ranging from 2 - 2000 M^r, are usually only slightly 
supercritical with < 5 ( Falgarone et al. 2008 ; Girart et al. 



2009 



Beuther et al.|2010 ) . Together with our numerical re- 



sults this suggest that typically there should be no well- 
collimated jets in the very early stages of high-mass star 
formation but rather sphere-like, slowly expanding outfiow 
structures. As discussed in the paper, the question of how 
well collimated jets are generated breaks down to the prob- 
lem of how Keplerian discs can be formed in highly magne- 
tised cores. For this to happen, effects in the later evolution 
of the system are required that reduce the efiaciency of mag- 
netic braking. In particular, the initial, slowly expanding 
outfiows as well as non-ideal MHD effects could accomplish, 
this resulting in a successive growth of discs and the devel- 
opment of well collimated jets over time. 
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